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We study Kosterlitz-Thouless (KT) transitions of the Larkin-Ovchinnikov (LO) phase for a two- 
dimensional system composed of coupled one-dimensional tubes of fermions. The low energy excita- 
tions involve the oscillation of the LO stripe and the fluctuation of the phase, which can be described 
by an effective theory composed of two anisotropic XY models. We compute from a microscopic 
model the coefficients of the XY models from which the KT transition temperatures are determined. 
We found the T KT oc t± for small intertube tunneling tj_ and over a wide range of parameters the 
fractional defects costs least energy to proliferate. As t± increases the system undergoes a first-order 
transition to the normal phase at zero temperature. Our method is efficient and can be used to 
determine the Goldstone excitations of any stripe order involving charge or spin degrees of freedom. 
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The role of topological excitations of striped superconducting states has been intensively studied [lH5[ since at 
finite temperature the proliferation of those defects can lead to possible exotic phases, such as the charge 4 superfluid 
5, 6]. A typical striped superconducting state is the Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) @, H[ state which is 
believed to exist in heavy-fermion superconductor CeCoIns [j| [l(| and has been recently proposed to occur in the 
system invol ving p-orbital bands [ill Il2j . Since the FFLO order is more likely to occur in the quasi-one-dimensional 
(ID) system [ijj, the cold atom system with two imbalanced species of atoms confined in a lattice array of ID tubes 
formed by coherent laser beams [14j seems more promising to display the direct evidence. Since the intertube coupling 
can be tuned relatively with ease in the cold atom system by controlling the intensity of trapping lasers, it is suitable 
to study the dimensional crossover phenomena [15H18l |. 

Numerous exotic phases have been predicted from effective field theories Q , but the phase diagram of these exotic 
phases is not established for cold atom experiments yet. In cold atom experiments, the microscopic parameters (like 
interaction strength) are tunable and measurable, and this motivates our detailed study of the Kosterlitz-Thouless 
(KT) transitions of the Larkin-Ovchinnikov (LO) phase starting from a microscopic model. In this work we study a 
quasi-lD two-dimensional (2D) spinful fcrmionic system composed of coupled ID tubes as illustrated in Fig. [T^ where 
at zero temperature the LO order is the ground state. We determine the Kosterlitz-Thouless temperature of LO 
phase (the FFLO regime in Fig. Q~|d 13]) as a function of intertube coupling t± from a microscopic model. We found 
that over a wide range of parameters the fractional defect dominates the phase transition and transition temperature 
is linear in tj_ for small t± (Fig. [5^). At zero temperature the transition from LO to normal phase, driven by the 
disappearance of the Fermi surface nesting upon increasing t±, is of first order (Fig. HJd). 

Before introducing the microscopic model, we first discuss the Goldstone modes of LO phase from the symmetry 
point of view. The LO phase is characterized by an order parameter of stripe configuration 
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FIG. 1: (a) Configuration of the system: arrays of ID tubes with intertube distance b and intertube tunneling t±. (b) A 
schematic plot of quasi-lD phase diagram as a function of fi and h. Vac: vacuum state (no particle); P-FL: partially polarized 
Fermi liquid; F-FL: fully polarized Fermi liquid. Our study here focuses on the FFLO regime, (c) Possible phases as a function 
of temperature. N:normal Fermi liquid; CDW: charge density wave; SF 4 : charge 4 superfluid. 
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FIG. 2: (Color online) (a) Phase diagram for /i = 2, h = 1.145, Q = 1.2, Ao = 0.22. Curves with ligands are computed whereas 
two dashed lines embracing the phase coexisting region are schematic. The calculated transition temperature T within the 
phase coexisting region is not well defined and its plot is switched to the dotted line with circles, (b) Energy as a function of 
gap amplitude A for y, = 2, h = 1.145, Q = 1.2, t± = — 0.2. The energy minimum occurs at A = 0.22 for t± < 0.15. 



where Ao is the amplitude of the order parameter while f(z) describes the stripe satisfying f(z + L/2) = — f(z) and 
f(z) — f(—z) in our coordinate choice. The LO wave vector Q is defined as Q = 2n/L with L the period of the stripe. 
Since the LO phase breaks both translational and U{1) symmetries, it has two branches of Goldstone modes - the 
oscillation of the stripe, and the phase fluctuation of the amplitude. Under these fluctuations, the order parameter 
becomes 

A(x,z)^A e^ x ^f(z + u(x 7 z)) (2) 

where u(x,z), 9{x,z) are generalized elastic fields [19] to describe the Goldstone modes. Physically u represents the 
small oscillation of the stripe LO order whereas 9 the phase fluctuation of the amplitude. In the quasi-lD system, 
x and z directions are not equivalent. Therefore to the quadratic order the total free energy in terms of generalized 
elastic fields are described by two anisotropic XY models [a, |20j 



AF = / dxdz 



(3) 



where f x — d x f. In our notation, u and 8/Q have the dimension of length, their first derivatives are dimensionless, and 
coefficients A, £?, C, D have the dimension of energy. For results presented in this paper, we take f{z) = cosQz which 
is (e l Q z +e~ l ® z )/2. In this case, these two Goldstone modes correspond to phase fluctuations of two Fulde-Ferrell (FF) 
order A(x,z) — A (e lQ ^ z+u ' + e ~*Q( z + u ))/2. For the FF order, A(x, z) oc e lQz which does not break translation 
symmetry. When identifying u + = u + 9/Q, u~ = u — 9/Q, A(x, z) — A cos(z + u(x, z)) e l6 ( x ' z )^ consistent with the 
functional form in Eq. ©. To be general, we shall not specify the form of f(z) unless necessary. 

In 2D, each elastic field is associated with one topological defect. For u the defect is the (edge) dislocation satisfying 
§ Vm • dl = Lrid', for 9 the defect is the vortex satisfying <f X79 ■ dl = 2itn v with n^, n v integers. There is another 
topological defect referred to as a half-vortex half-dislocation (HH) where {rid,n v ) — (±1/2, ±1/2), which originates 
from the Z2 symmetry of the order parameter (HQ - when circulating around an HH defect, each of the half vortex and 
the half dislocation introduces a minus sign leaving the order parameter unchanged. The proliferation of topological 
defects leads to Kosterlitz-Thouless (KT) transition. The KT transition temperature T u = ^ V AB for dislocations, 
T e = %VCD for vortices, and T hh = f (VAB + VCD) = (T u + T e )/A for HH [5]. The last temperature cannot be 
highest. When increasing the temperature, there are three distinct possibilities as illustrated in Fig. [TJ;: (I) T hh is 
the lowest, (II) T e the lowest, and (III) T u the lowest For (I) there is only one transition from LO to normal state 
at T hh . For (II), the LO phase first becomes a charge density wave (CDW) state at T 9 and then normal at T hh . For 
(III), the LO phase first becomes a charge 4 superfluid at T u and then normal at T hh . 

The microscopic model we use is a one band model with attractive contact interaction in a quasi-lD system [l3l j. 

H = V^r ct c r 

k,cr 

+ 9iD^2 dzc l x ,t( z ) c Li( z ) C4 -^( z ) Cl -t( z ) ( 4 ) 



2, ,2 



ft k~ 

with £7; = -7^- — 2t_iCOs(bk x ) — fj, + h(— l) a , where k z is unbound, \bk x \ < it with b the intertube distance, and 
(— l)t = 1, (—1)^ = —1. Following Ref [H, [2l|, we measure all lengths in the ID scattering length a ±D — — 
and all energies by the ID bound energy e_B = ■ The dimensionless parameters in this model are tj_/eB, fJ-f^B, 
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h/eB, and b/a\D- The attractive interaction implies negative q\d and thus positive anj. The relation between a\D 
and a 3 o is 123 a id — ( 1,4 ^ )3 ) with aj_ = W —7— . The typical is of the order lOOnm and a 3 £> can be 

\ "3D v 2 / y <nu/j_ 

controlled by the Fashbach resonance. Taking ci\d = lOOnm, m =6/(6 x 10 23 )g ( 6 Li) [l4j]) the bound state energy is 
e B ~ 1.5 x 1(T 6 K. 

To obtain the effective theory at given parameters, we first solve the Hamiltonian by the variational method where 
the order parameter is assumed to be sinusoidal A i:c {z) = A{z) = A cos(Qz) = 9iD{ c i x ,l(z)ci K ,f(z))- f(z) in Eq. @ 
is chosen to be cos Qz. The mean field (Bogoliubov de Gennes) Hamiltonian is 



where 



H (A„, C )=(% (6) 

with the block B = &k x ,~k x {$kz+Q ,-k z + 3k z -Q,-k z )' The Ao and Q are determined by minimizing the free energy 
with respect to Q and Ao given by 

F[A] = -T5>g[l + e^] + + L * N *^1 (7) 

where e„ are a/1 eigenvalues of the matrix if (Ao, Q). In the T — > limit, — Tj^ n log[l + e _e "/ T ] = e n G(— e„). 

The energy cost for given strain configurations u(x, z), 9(x, z) of elastic fields is computed by AF(b) — F[A^e l9 f(z + 
u)} — F[Aof(z)}, which to the lowest order of u and 9 in quasi-lD reduces to the form 



AF(b) = bJ2j 



dz 



2 

C{b) o2 , D(b) a2 



A{b \Qu x f + ^(Qu z f 



(8) 



where in the quasi-lD system u x (xi) = u ( x i+ b )~ u ( x i) ; 9 x {xi) = £i2»ijh^fe) i Replacing 6^ — > / dx, Eq. © is the 
same as Eq. ([3]). 

We now show that the KT transition temperature is independent of the intertube distance b when t± is fixed. For 
simplicity we only consider the stripe oscillation field u. The same argument applies to the phase field 9. When 
b — > ab, the total free energy cost due to compression or stretching along z is unchanged (since the energy depends 
only on t± which is fixed) but the energy density changes. Consequently AF z (b) — AF z (ab) implies 

bN x L z ^l(Qu z ) 2 = a bN x L z ^^{Qu z f (9) 

leading to B(ab) = B(b)/a. The free energy caused by u x depends only on intertube coupling t± and therefore the 
elastic field difference between two adjacent tubes u(x + b) — u{x). When b — > ab, as long as u(x + b) — u(x) = 
u(x + ab) — u(x) (same t±) the total free energy remains unchanged which leads to 

= ^l Lziab) srQ 2 ^(* + <^ (10) 

which leads to A(ab) = aA(b). Their product A(b)B(b) is independent of b, so is the transition temperature T KT oc 
\/ AB. Since our main interest is the transition temperature, we take b = aiD = 1. 
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To obtain coefficients A, B, C, D in Eq. ([3]), we take the following approach. Take B as an example, we choose 
u(x,z) = u z z, 9(x,z) = 0, compute AF(b) for several u z , and fit SF(b,u z ) = AF(b; u z ) / (bN x L z ) = ^rp-(Qu z ) 2 as 
expressed in Eq. ©. The same procedure apply to A, C, D. There is another approach to obtain these coefficients 
involving Green's function [23] which requires computing the inverse of a matrix and is very time-consuming. Our 
approach instead only involves the computation of eigenvalues [24j which allows us to include more fc-points. 

Some key steps of computing 6F(b, u z ) are are summarized here. To obtain the coefficients C and D requires 
computing the eigenvalues of 



H(0) 



(h - fx) - 
A /(z)e- 



A f(z)e ie 
-(h - £t) +h 



(11) 



h k 

where ho is diagonalizcd in momentum space as — 2t± cos bk x . By performing a local gauge transformation 
c(z) —> c(z)e~ ie l 2 , H(9) in new coordinate becomes Hi(9) 0, [25| which is 



Hl (6) = ( , Ao W 

1 A f(z) -/la,; 



(12) 



€-n+*/2,i respectively 



where ho t -f and ho,], are diagonalized in momentum space as £ + (fc) = ££ + #/2 f anc l £ (~ ^) 

with v — V0/2. The eigenvalues of H and H\ are identical, but this transformation automatically obtains the 
derivative of 8, i.e. 8 X and 8 Z , in the diagonal blocks. We use fx = 2, h = 1.145, t± = 0.1 as an example. Minimizing 
the energy functional with respect to A and Q leads to Q — 1.2, A = 0.22. Fig. [3] shows 8F(9 Z ) for 9(x,z) = 9 z z 
and SF(9 X ) for 9(x, z) = 9 x x from which the quadratic fit leads to C = 0.00168 and D = 0.23. 





FIG. 3: (Color online) Energy density as a function of 9 Z and 9 X for fi — 2, h = 1.145, Q = 1.2, Ao = 0.22.. 



To obtain the coefficients A and B requires computing the eigenvalues of 

A useful trick is to do the calculation in a new coordinate whose order parameter is exactly Ao/(2)[23j. We stress here 
the Jacobian arising from the coordinate transform has to be taken into account because it is the free energy, not the 
free energy density, which is invariant under coordinate transformation. Again we use parameter /i = 2, h — 1.145, 
t± = 0.1 as an example. Fig. |4] shows SF(Qu z ) for u(x, z) = u z z and 5F(Qu x ) for 9(x. z) — u x x from which we can 
fit A = 0.00168 and B = 0.234. Note that f(z) is taken to be cosQz for the results presented here. However we 
emphasize that the coefficients A, B, C, D can be obtained for any given order parameters. 




FIG. 4: (Color online) Energy density as a function of u z and u x for /i = 2, h — 1.145, Q = 1.2, Ao = 0.22. 
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We compute the coefficients for severai t± and determine aii three KT transition temperatures. Our main result is 
shown in Fig. [5^ where the phase diagram as a function of T and ij_ is plotted for a representative set of parameters 
fi = 2, h = 1.145, Q = 1.2, Ao = 0.22. When t± = 0, there is no intertube coupling and no correlation along x 
leading to zero T KT . As the system goes from pure ID to quasi- ID, T KT (t±) oc t±. More specifically, we found the 
coefficients associated with z derivative, i.e. B and D, depend very weakly on t±_ whereas those with x derivative, 
i.e. A and C, depend quadratically on t±. This explains the linear t± dependence of T KT (oc s/AB, y/CD). We 
found T e and T u are very close because the "cosine" ansatz is very close to two decoupled FF order with opposite 
wave vector for small Ao and it is the coupling between u + and u~ (the fluctuations of two FF orders) which lifts the 
degeneracy of u and 9 fields. In this case the proliferation of half-vortex half-dislocation costs the least energy and 
the system only undergoes one transition from LO to normal at T hh when raising the temperature, as shown in the 
case I of Fig. QJ. The transition temperature is of the order of 0.03 cb which is roughly 8 x 10 _8 K for the system 
of 6 Li with ain = lOOnm. One also notes that the obtained T KT (~ 0.03) is an order of magnitude smaller than the 
mean field gap Ao(= 0.22), so the coefficients computed at T = are almost identical (less than 1% difference) to 
those computed at T ~ T KT . 

At T = 0, our simulation suggests the quantum phase transition from LO to normal phases upon increasing t±_ is 
of first order. Fig. ^Bp shows the energy as a function of A for /j, = 2, h = 1.145, Q = 1.2, t± £ (0,0.2) where the 
minimum determines the value of Ao. We found that as t± increases Ao stays around 0.22 and when t±_ > 0.16 Ao 
becomes zero. Around t± = 0.16, the E(A) is essentially flat with several shallow minima. We note that the FFLO 
to normal transition as a function of temperature (fixed fi, and h) 0, [13] or h (fixed [i, T) (26| is also of first order. 
At finite temperature around the critical t±, the system is in the phase coexisting region. 

We have assumed a sinusoidal order parameter in the current calculation. However near BCS/LO transition (Fig. 
Eb) [HLHl]; the order parameter behaves more domain- wall like (2(| [53] than sinusoidal. Therefore the ansatz with 
sinusoidal order parameter does not capture all physics. Close to the BCS/LO transition, we expect that the stripe 
fluctuation should be stronger than the phase fluctuation (A < C, B < D) and a two-stage transition with charge 
4 superfluid shown as the case III in Fig. QF. can happen. In a cold atom trap where the chemical potential is a 
position-dependent, the interface between phases shown in Fig. [TJ; is unavoidable and worth investigating. 

In summary, we have computed from a microscopic model the effective theories of Goldstone modes of the LO 
order for a quasi-lD fcrmionic system from which the Kosterlitz-Thouless transition temperatures are determined. 
The transition temperatures are found to depend linearly on the intertube coupling t±. However, the method applied 
here ne glec ts the quantum fluctuation along the ID tubes which can modify this linear t± dependence, especially at 
t_L — > [23]. As t± increases, the system goes to a phase coexisting regime sandwiched by the LO and normal phases. 
Our approach can generally determine the Goldstone excitations of any stripe order involving charge or spin from a 
microscopic model which should be useful for comparison between theories and experiments. 
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No. PHY05-51164. 
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